Possible chiral structure in the pseudogap phase of cuprates 
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We propose a novel chiral order parameter to explain the unusual polar Kerr effect in underdoped 
cuprates. It is based on the loop-current model by Varma, which is characterized by the in-plane 
anapole moment N and exhibits the magnetoelectric effect. We propose a helical structure where 
the vector JV^™- 1 in the layer n is twisted by the angle tt/2 relative to iV"'™ -1 ', thus breaking inversion 
symmetry. We show that coupling between magnetoelectric terms in the neighboring layers for this 
structure produces optical gyrotropy, which results in circular dichroism and the polar Kerr effect. 

PACS numbers: 74.72.Gh, 78.20.Jq, 78.20.Ek 
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Introduction. The nature of the pseudogap phase in 
underdoped cuprate superconductors has been a long- 
standing problem []]]. A series of optical measurements 
0-01 revealed gyrotropy in this state. It was observed 
that the polarizations of incident and reflected light dif- 
fer by a small angle 9k, called the polar Kerr angle. Ini- 
tially, these experiments were interpreted as the evidence 
for spontaneous time-reversal symmetry breaking. The- 
oretical models derived optical gyrotropy from the 
anomalous Hall effect. In these scenarios, the order pa- 
rameter is equivalent to an intrinsic magnetic field per- 
pendicular to the layers, which permeates the system and 
points inward and outward at the opposite surfaces of a 
crystal. Therefore, the Kerr angle should have opposite 
signs at the opposite surfaces of the crystal. 

However, recent reports 1(J II | found that the Kerr 



angle has the same sign at the opposite surfaces of a sam- 
ple. Therefore, the observed gyrotropy is not consistent 
with the time-reversal-symmetry breaking due to a mag- 
netic order and should be interpreted as the evidence 
for natural optical activity due to inversion symmetry 
breaking [12j. Systems with helical structures, such as 
cholesteric liquid crystals and some organic molecules, 
typically exhibit optical gyrotropy and the polar Kerr ef- 
fect. It is important that the sign of the Kerr angle in this 
case is the same at the opposite surfaces of the system, in 
contrast to the gyrotropy produced by a magnetic order 
(see Ch. 11 and 12 in Ref. " 



13) 



Recent papers [llj, [13J proposed various theoretical 



scenarios for appearance of a chiral order parameter in 
cuprates. Three possible orders were studied phenomeno- 
logically in Ref. [llj: 3-plane rotation of distorted Fermi 
circles, 4-plane rotation of a structure with an in-plane 
ferroelectric moment, and 4-plane rotation of an in-plane 
density wave with period 3. A formula for gyrotropy of 
chiral metals was derived in Ref. (Hf in terms of the Berry 
curvature in the momentum space. However, Ref. [14| 
questioned applicability of the implied assumption of the 
coherent electron motion between the layers in Ref. [l3[ . 

Helical order of the anapole moments. Here we pro- 
pose a novel chiral state, which could account for the 



polar Kerr effect in cuprates without requiring coherent 
electron motion between the layers. The starting point 
for our construction is the model of persistent loop cur- 
rents by Simon and Varma (l5j shown in Fig. Ufa). The 
configuration of persistent currents is such that the to- 
tal magnetic flux through the lattice unit cell vanishes. 
The anomalous Hall effect is zero [l6| , and this state does 
not exhibit magnetic gyrotropy. The loop current order 
[13] is characterized by the anapole moment N 18, 19j . 
shown by the green arrow in Fig. HJa) and defined as 



N = 



J d 2 r[m(r) x r] = ^- J d 2 rr 2 j(r). (1) 



Here m{r) and j(r) are the microscopic densities of the 
magnetic moment and electric current, and the integral 
is taken over the unit cell. The square symmetry of the 
lattice allows four possible orientations [20j for the vector 
AT, which can be obtained by progressive tt/2 rotations 
of the configuration shown in Fig. [lja). 

We propose a novel chiral state where the anapole mo- 
ments iVW in consecutive CuC>2 layers labeled by the 
integer number n are rotated by tt/2, so that they trace 
out a helix, as shown in Fig. [TJb). This arrangement is 
somewhat similar to cholesteric liquid crystals [21] . This 
spiral structure breaks three-dimensional inversion sym- 
metry and can be qualitatively visualized as follows. The 
in-plane triangular loop currents in Fig.QJa) produce per- 
pendicular magnetic fields of the opposite signs shown by 
the circled red dot and blue cross at the centers of the 
loops. When the anapole moments N^ n > are arranged in 
the spiral structure in Fig. QJb), the red and blue mag- 
netic field lines, propagating from one layer to another, 
form a double helix due to twisting of N^ n K 

The chiral structure is characterized by a pseudoscalar 
order parameter a changing sign upon inversion 



a 



\N( n ) x JV < " +1) ' 



(2) 



where z is the unit vector along the z axis. If the anapole 
vector 7V( n ) is static, Eq. ([2|) does not need the brackets 
for averaging. However, if the vector N^- 71 ' fluctuates, the 
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FIG. 1: (Color online) (a) Loop-current order in a CuC>2 plane 
[l5l ]. Black arrows show directions of microscopic persistent 
currents between copper and oxygen atoms. Green arrow 
shows the anapole moment N . (b) Chiral order constructed 
on a series of parallel Cu02 planes. The vector N rotates by 
the angle 7r/2 from one layer to another, and the period of 
the structure is fourfold. The blue and red curves are the two 
magnetic field lines that intertwine in a double helix. 

brackets in Eq. J2J represent thermodynamic and, pos- 
sibly, quantum |22j averaging. The chiral order param- 
eter a is a local correlation function of the orientations 
of 7V(™) in the neighboring layers and does not require 
long-range order in N^ n ' . Spontaneous chiral symmetry 
breaking is known for other systems 0, [24| ■ The config- 
uration with 7r/2 rotations in Fig.QTb) maximizes a for a 
given absolute value of N. Notice that Eq. ([2]) is similar 
to the Dzyaloshinskii-Moriya interaction for spins. 

In electrodynamics of media 12], natural optical ac- 
tivity arises when inversion symmetry is broken and the 
expansion of the dielectric tensor e^ v {uj,k) in powers of 
the wave vector k has a non-vanishing first-order term 



■ lv (u), k) = e MV (u, 0) + ij(u>)e^ z k z . 



(3) 



Here e M „A is the antisymmetric tensor, and k z is the wave 
vector of an electromagnetic wave propagating along the 
z axis. The second term in Eq. ^ represents a non-local 
effect along the z axis and is responsible for gyrotropic 
properties of the medium. The polar Kerr angle 9k is 
determined by the following formula [lH to the first order 
in the coefficient 7 in Eq. <j3j> 

7(w) 



tan#K-(a;) = —Xm 



1 - e(w) 



(4) 



It is clear that a non-zero Kerr angle requires an imagi- 
nary part, i.e. dissipation, either in s(u) or y(ui). In the 
rest of the paper, we derive the second term in Eq. Q 
for the spiral structure in Fig. [TJb). 
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FIG. 2: (a) Left: Feynman diagram for the effective action 
of electromagnetic fields (wavy lines), obtained by integrating 
out the electron field (solid lines with arrows). Right: The 
magnetoelectric term in the effective action, Eq. where 
the double line represents the anapole moment N . (b) Cou- 
pling between the magnetoelectric terms at the neighboring 
layers produces the effective action for the electric field in 
Eq. ©. The dashed wavy line represents the magnetic field 
propagator, and the double solid line represents the interlayer 
correlator of the anapole moments in Eq. ([2]). 



Magnetoelectric effect in a single C11O2 plane. First 
let us consider a single Cu02 plane with loop currents in 
Fig. [IJa). Integration of the electron field, schematically 
shown in the left diagram in Fig. Ufa), gives an effective 
action for the electromagnetic field with a magnetoelec- 
tric term [l2|, . By symmetry, it has the form 17 1 



Sme = / dujd 2 rp{uj)B z {~u:) [N x E{u)] g . (5) 



Here E = (E x ,E y ) is the in-plane electric field, B z is 
the out-of-plane magnetic field, and N = (N x , N y ) is the 
in-plane anapole moment. Given Eq. ([1]), the anapole 
moment N is the time-reversal-odd and parity-odd vec- 
tor, so Eq. ([5]) has the correct symmetry structure. It is 
represented graphically by the right diagram in Fig.f2Ja). 
Equation ([S]) is written in the frequency representation 
for the electromagnetic fields, whereas TV is taken to be 
static, i.e. having zero frequency, and P(bS) is a frequency- 
dependent coefficient. The effective action in Eq. ([5]) is 
written in the continuous, long- wavelength limit by av- 
eraging the electromagnetic fields over distances longer 
than the unit cell of the lattice. 

By taking a variation of Eq. ([5]), we find that an in- 
plane electric field induces an out-of-plane magnetization 



M z (u) 



SS- 



ME 



5B z (-lo) 



= [N x E{u>)] z . 



(6) 



Physical interpretation is clear by symmetry in Fig. [TJa) . 
An in-plane electric field E _L N promotes electron trans- 
fer from one triangular loop to another, thus breaking 
symmetry and making one loop current stronger, which 
results in the net perpendicular magnetization. 

Similarly, an out-of-plane magnetic field induces an in- 



3 



plane electric polarization 
SSme 



P(u) = 



SE(-uj) 



= 0{u) [NxB z (u)} 



(7) 



The perpendicular magnetic field B z lowers the energy 
for one loop current and increases for another in Fig.[lja), 
which results in electron transfer between the loops and 
the in-plane electric polarization P _L TV. 

Interlayer coupling and the gyrotropic term. Now let 
us consider many parallel layers labeled by the integer 
number n. In this case, we expect a magnetic coupling 
between magnetizations at the neighboring layers 



Smm = 



duj d 



E 



AM^ n \uj) M^ n+1 \—oj). (8) 



This term should be derived microscopically from the 
distribution of electric currents inside the unit cell com- 
parable to the interlayer distance d. Here we simply 
write Eq. flSJ by symmetry for the long- wavelength fields 
M z with a phenomenological coefficient A. We assume 
that the interlayer coupling between magnetizations is 
instantaneous, i.e. A is frequency-independent. This as- 
sumption is valid when the interlayer distance d is much 
smaller than the wavelength of light: d -C A = 27r/fc 2 . 

Substituting Eq. (JB]) into Eq. ([5]) , we obtain an effective 
action for the electric fields in the multilayer system 

See = - f du:d 2 r\l3(Lo)\ 2 Y,k (9) 

J n 

[JVW x E^(lu)] z [N in+1 ^xE^ +1 \~uj)] z , 

where we used the standard relation j3(—u>) — for 
a linear response function. Figure [HL) illustrates this 
calculation diagrammatically. By coupling the magneto- 
electric vertices shown in Fig. [U[a) and integrating out 
the magnetic field propagator shown by the dashed wavy 
line, we obtain the effective action for the electric field in 
Eq. (|9|). The double solid line represents the interlayer 
correlator of the anapole moments in Eq. ([2]). 

Let us choose the x and y axes along the crystallo- 
graphic direction a + b and b — a in Fig. [TJ so that the 
vectors = — A/"( n+2 ) are along x for odd n and y 

for even n. Expanding the vector products in Eq. we 
find two terms in the sum, for odd and even n. Changing 
the variable n — > n — 1 in the latter sum, we find 

See- / dLud 2 r\P(u)\ 2 AN^N^+V 



4 n) M [E^H-^-Elp-Vi-u,)] 



(10) 



Using Eq. ^ , taking the continuous limit z = nd, where 



(n+l) 



= 2d{dE x /dz) and 2d£ nodd = f dz, 



E. 

and integrating by parts in z, we get 



Ske = ~ / dwd*r\p{Lj)fz 



E(uj) x 



8E(-oj) 
dz 



or equivalently 

See = ^ Jdud 3 r\p{uj)\ 2 E(Lj)-[V z xE(-Lo)}. (12) 
Comparing Eq. (1111) with the standard expression 

S =^J dujd 3 k£» v (u,k)E l ,(u,k)E u (-Lj,-k), (13) 

we obtain the coefficient 7 in the second term in Eq. ([3]) 

7 (w) = 4ttqA|/3(w)| 2 . (14) 

Equation (|14[) shows that the gyrotropic coefficient 7(0;) 
is determined by the chiral order parameter a, the in- 
terlayer magnetic coupling A, and the magnetoelectric 
coefficient /3(w). The sign of 7 depends on the sign of a. 

The above derivation was presented for equally spaced 
Cu02 layers. However, many cuprates have the bilayer 
structure, where the interlayer distances alternate be- 
tween d =F Ad. In this case, the interlayer coupling co- 
efficient in Eq. ([8]) is A 1 ™) = A ± £ for even and odd n. 
As a result, we find an additional term, which is similar 



to Eq. pli]). but with A -( and E x ' 



(n+l) 



■Ei n - 1) -> 



Er^'+E^ L > = 2E X . This term contributes an off- 
diagonal symmetric term to the dielectric tensor e xy = 
e yx = — iTraC,\f3(uj)\ 2 /d, which becomes diagonal in the 
crystallographic basis of a and b 



-e fcb = 47raCI/3H| 2 M 



(15) 



(11) 



Thus, we find that the helical structure in the presence of 
bilayers produces nematicity, i.e. anisotropy between the 
crystallographic directions a and b. This is clear by sym- 
metry in Fig.JlJb), where the pairs of layers (n, n+l) have 
the preferred direction +N < ^ n+1 " > along a. Equations 
(fT4"|) and (fT5|) generate circular and linear dichroism. If 
A ~ £, the linear dichroism is much stronger than the 
circular one because cZ <gc A = 2n/k z . 

Figure [3] illustrates these results graphically. The elec- 
tric field component fiW !_ generates a magneti- 

(n) 

zation M z at the bottom layer via the magnetoelectric 
effect in Eq. ([6]). The magnetization Mz produces a 
magnetic field B z n+1 ^ in the same direction at the top 
layer via Eq. ©. This magnetic field induces an elec- 
tric polarization p(™ +1 ) J_ iV(" +1 ) at the top layer via 
the magnetoelectric effect in Eq. ([7]). Taking into ac- 
count the third layer n + 2 (not shown), we find that 
p(n+i) ^ £ x [ C2 _e(»+2) _ Cl #(«■)] with some coefficients 

ci and C2- For equally spaced layers with c\ = C2, we get 
Eq. (fTT|) . For bilayers with c\ ^ C2, we get the additional 
nematic term in Eq. (|15|) . 

Experimental relevance. Let us discuss possible ex- 
perimental evidence for the proposed chiral order in 
cuprates. Neutron scattering measurements [27, 28 1 pro- 
vide support for the loop currents shown in Fig. [lja). 
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FIG. 3: (Color online) Schematic illustration of the rotation 
of electric polarization for the chiral state shown in Fig. [lib). 
The electric field at the bottom layer induces the mag- 

netization M^ n \ which couples to the top layer and induces 
the electric polarization p( n+1 > _|_ because of the twist 

in the anapole moments iV (n) ± jV (n+1) . 



However, the NMR and NQR experiments [29] find no 
evidence for the local magnetic fields predicted by this 
model. So, the experimental situation remains contro- 
versial. A survey of experimental evidence supporting 
loop currents is presented in Ref. f3Tjj |. 

Although the loop current order in Fig. [JJa) breaks ro- 
tational symmetry in the plane as specified by the vector 
TV, the neutron scattering measurements [23, [28| always 
observe the full rotational symmetry. This may be due 
to domains with different TV, but the spiral order shown 
in Fig. UJb) also provides a natural explanation. While 
the system has the tendency to break rotational symme- 
try in each CuC>2 layer, we argue that it tries to restore 
macroscopic symmetry by orienting the vectors TV or- 
thogonally in the neighboring layers, which is consistent 
with the spiral structure. 

Moreover, Ref. [27]] observed that the microscopic mag- 
netic field from by the loop currents is not perpendicular 
to the layers, as expected from Fig. [JJa), but has an in- 
plane component. This effect can be explained by the 
spiral order in Fig. UJb) [3l|. Since the magnetic field 
lines are twisted in a double-helix structure, they are 
naturally tilted with an in-plane component. In the pres- 
ence of TV, the energy of the system contains the term 
TV • [V x B] = B ■ [V x TV] Since it is linear in 

B, whereas magnetic energy is proportional to B 2 , the 
system develops an equilibrium in-plane magnetic field 
B oc [V x TV] parallel to TV in Fig. [Tp, in qualitative 
agreement with Ref. [27]. Moreover, the total energy de- 
creases as — |V x TV| 2 , which favors the spiral structure. 



Finally, the recent X-ray measurements (32H34j found 
doubling of the unit cell in YBCO in the z direction [3o| . 
Given the bilayer structure of YBCO, the new unit cell 
contains four Cu02 layers. The fourfold period is consis- 
tent with the spiral order shown in Fig. [Tib). 



For bilayer materials, the spiral structure in Fig. [Tfb) 
naturally produces nematicity, where the oxygen atoms 
to the left and right of the copper atom in Fig. QJa) are 
not equivalent to the oxygen atoms above and below. 
This nematic symmetry is in qualitative agreement with 
the pattern observed in the scanning tunneling measure- 
ments |36| . although the same pattern was observed ex- 
perimentally in bilayer and single-layer cuprates. 

Conclusions. We propose a fourfold chiral state for 
cuprates obtained by twisting Varma's loop current or- 
der by 7r/2 in consecutive Cu02 layers. We show that 
this state exhibits natural optical activity and derive the 
gyrotropic coefficient. It can account for the polar Kerr 
effect in cuprates [2-5] without invoking magnetic gy- 



rotropy [10 . For bilayer compounds, we also find 



nematicity and linear dichroism. Our model is based on 
magnetic coupling between the C11O2 layers and does not 
require coherent electron tunneling between the layers 
and long-range order in the chiral structure. Other mod- 
els for the polar Kerr effect in cuprates invoked magne- 
toelectric effects 13711 and coupling between loop currents 
with different TV [lfjj, but considered only a single layer, 
rather than the spiral multilayer structure. 
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